EXPANDER – an integrative program suite for microarray data analysis

BACKGROUND: Gene expression microarrays are a prominent experimental tool in functional genomics which has opened the opportunity for gaining global, systems-level understanding of transcriptional networks. Experiments that apply this technology typically generate overwhelming volumes of data, unprecedented in biological research. Therefore the task of mining meaningful biological knowledge out of the raw data is a major challenge in bioinformatics. Of special need are integrative packages that provide biologist users with advanced but yet easy to use, set of algorithms, together covering the whole range of steps in microarray data analysis. RESULTS: Here we present the EXPANDER 2.0 (EXPression ANalyzer and DisplayER) software package. EXPANDER 2.0 is an integrative package for the analysis of gene expression data, designed as a 'one-stop shop' tool that implements various data analysis algorithms ranging from the initial steps of normalization and filtering, through clustering and biclustering, to high-level functional enrichment analysis that points to biological processes that are active in the examined conditions, and to promoter cis-regulatory elements analysis that elucidates transcription factors that control the observed transcriptional response. EXPANDER is available with pre-compiled functional Gene Ontology (GO) and promoter sequence-derived data files for yeast, worm, fly, rat, mouse and human, supporting high-level analysis applied to data obtained from these six organisms. CONCLUSION: EXPANDER integrated capabilities and its built-in support of multiple organisms make it a very powerful tool for analysis of microarray data. The package is freely available for academic users a

Predicting risk for Alcohol Use Disorder using longitudinal data with multimodal biomarkers and family history: a machine learning study.

Predictive models have succeeded in distinguishing between individuals with Alcohol use Disorder (AUD) and controls. However, predictive models identifying who is prone to develop AUD and the biomarkers indicating a predisposition to AUD are still unclear. Our sample (n = 656) included offspring and non-offspring of European American (EA) and African American (AA) ancestry from the Collaborative Study of the Genetics of Alcoholism (COGA) who were recruited as early as age 12 and were unaffected at first assessment and reassessed years later as AUD (DSM-5) (n = 328) or unaffected (n = 328). Machine learning analysis was performed for 220 EEG measures, 149 alcohol-related single nucleotide polymorphisms (SNPs) from a recent large Genome-wide Association Study (GWAS) of alcohol use/misuse and two family history (mother DSM-5 AUD and father DSM-5 AUD) features using supervised, Linear Support Vector Machine (SVM) classifier to test which features assessed before developing AUD predict those who go on to develop AUD. Age, gender, and ancestry stratified analyses were performed. Results indicate significant and higher accuracy rates for the AA compared with the EA prediction models and a higher model accuracy trend among females compared with males for both ancestries. Combined EEG and SNP features model outperformed models based on only EEG features or only SNP features for both EA and AA samples. This multidimensional superiority was confirmed in a follow-up analysis in the AA age groups (12-15, 16-19, 20-30) and EA age group (16-19). In both ancestry samples, the youngest age group achieved higher accuracy score than the two other older age groups. Maternal AUD increased the model's accuracy in both ancestries' samples. Several discriminative EEG measures and SNPs features were identified, including lower posterior gamma, higher slow wave connectivity (delta, theta, alpha), higher frontal gamma ratio, higher beta correlation in the parietal area, and 5 SNPs: rs4780836, rs2605140, rs11690265, rs692854, and rs13380649. Results highlight the significance of sampling uniformity followed by stratified (e.g., ancestry, gender, developmental period) analysis, and wider selection of features, to generate better prediction scores allowing a more accurate estimation of AUD development

Assessment of brain age in posttraumatic stress disorder: Findings from the ENIGMA PTSD and brain age working groups

Background: Posttraumatic stress disorder (PTSD) is associated with markers of accelerated aging. Estimates of brain age, compared to chronological age, may clarify the effects of PTSD on the brain and may inform treatment approaches targeting the neurobiology of aging in the context of PTSD. Method: Adult subjects (N = 2229; 56.2% male) aged 18–69 years (mean = 35.6, SD = 11.0) from 21 ENIGMA-PGC PTSD sites underwent T1-weighted brain structural magnetic resonance imaging, and PTSD assessment (PTSD+, n = 884). Previously trained voxel-wise (brainageR) and region-of-interest (BARACUS and PHOTON) machine learning pipelines were compared in a subset of control subjects (n = 386). Linear mixed effects models were conducted in the full sample (those with and without PTSD) to examine the effect of PTSD on brain predicted age difference (brain PAD; brain age − chronological age) controlling for chronological age, sex, and scan site. Results: BrainageR most accurately predicted brain age in a subset (n = 386) of controls (brainageR: ICC = 0.71, R = 0.72, MAE = 5.68; PHOTON: ICC = 0.61, R = 0.62, MAE = 6.37; BARACUS: ICC = 0.47, R = 0.64, MAE = 8.80). Using brainageR, a three-way interaction revealed that young males with PTSD exhibited higher brain PAD relative to male controls in young and old age groups; old males with PTSD exhibited lower brain PAD compared to male controls of all ages. Discussion: Differential impact of PTSD on brain PAD in younger versus older males may indicate a critical window when PTSD impacts brain aging, followed by age-related brain changes that are consonant with individuals without PTSD. Future longitudinal research is warranted to understand how PTSD impacts brain aging across the lifespan

Assessment of Brain Age in Posttraumatic Stress Disorder: Findings from the ENIGMA PTSD and Brain Age Working Groups

Background Posttraumatic stress disorder (PTSD) is associated with markers of accelerated aging. Estimates of brain age, compared to chronological age, may clarify the effects of PTSD on the brain and may inform treatment approaches targeting the neurobiology of aging in the context of PTSD. Method Adult subjects (N = 2229; 56.2% male) aged 18–69 years (mean = 35.6, SD = 11.0) from 21 ENIGMA-PGC PTSD sites underwent T1-weighted brain structural magnetic resonance imaging, and PTSD assessment (PTSD+, n = 884). Previously trained voxel-wise (brainageR) and region-of-interest (BARACUS and PHOTON) machine learning pipelines were compared in a subset of control subjects (n = 386). Linear mixed effects models were conducted in the full sample (those with and without PTSD) to examine the effect of PTSD on brain predicted age difference (brain PAD; brain age − chronological age) controlling for chronological age, sex, and scan site. Results BrainageR most accurately predicted brain age in a subset (n = 386) of controls (brainageR: ICC = 0.71, R = 0.72, MAE = 5.68; PHOTON: ICC = 0.61, R = 0.62, MAE = 6.37; BARACUS: ICC = 0.47, R = 0.64, MAE = 8.80). Using brainageR, a three-way interaction revealed that young males with PTSD exhibited higher brain PAD relative to male controls in young and old age groups; old males with PTSD exhibited lower brain PAD compared to male controls of all ages. Discussion Differential impact of PTSD on brain PAD in younger versus older males may indicate a critical window when PTSD impacts brain aging, followed by age-related brain changes that are consonant with individuals without PTSD. Future longitudinal research is warranted to understand how PTSD impacts brain aging across the lifespan

Remodeling of the Cortical Structural Connectome in Posttraumatic Stress Disorder:Results from the ENIGMA-PGC PTSD Consortium

BACKGROUND: Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA). METHODS: Neuroimaging and clinical data were aggregated from 29 research sites in >1,300 PTSD cases and >2,000 trauma-exposed controls (age 6.2-85.2 years) by the ENIGMA-PGC PTSD working group. Cortical regions in the network were rank-ordered by effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2 to 148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared to the mean SC of 5,000 randomly generated n-region networks. RESULTS: Patients with PTSD, relative to non-PTSD controls, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex and age modulated covariance differences of PTSD-related structural networks. CONCLUSIONS: Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The structural covariance networks that are perturbed in PTSD comport with converging evidence from resting state functional connectivity networks and networks impacted by inflammatory processes, and stress hormones in PTSD

A Comparison of Methods to Harmonize Cortical Thickness Measurements Across Scanners and Sites

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants’ demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2–81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3–85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ 2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ 2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ 2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects

Using fMRI connectivity to define a treatment-resistant form of post-traumatic stress disorder.

A mechanistic understanding of the pathology of psychiatric disorders has been hampered by extensive heterogeneity in biology, symptoms, and behavior within diagnostic categories that are defined subjectively. We investigated whether leveraging individual differences in information-processing impairments in patients with post-traumatic stress disorder (PTSD) could reveal phenotypes within the disorder. We found that a subgroup of patients with PTSD from two independent cohorts displayed both aberrant functional connectivity within the ventral attention network (VAN) as revealed by functional magnetic resonance imaging (fMRI) neuroimaging and impaired verbal memory on a word list learning task. This combined phenotype was not associated with differences in symptoms or comorbidities, but nonetheless could be used to predict a poor response to psychotherapy, the best-validated treatment for PTSD. Using concurrent focal noninvasive transcranial magnetic stimulation and electroencephalography, we then identified alterations in neural signal flow in the VAN that were evoked by direct stimulation of that network. These alterations were associated with individual differences in functional fMRI connectivity within the VAN. Our findings define specific neurobiological mechanisms in a subgroup of patients with PTSD that could contribute to the poor response to psychotherapy.PEV was supported by the Medical Research Council (grant no. MR/K020706/1) and is a Fellow of MQ: Transforming Mental Health (MQF17_24)

CLICK and EXPANDER: A System for Clustering and Visualizing Gene Expression Data

Motivation: Microarrays have become a central tool in biological research. Their applications range from functional annotation to tissue classification and genetic network inference. A key step in the analysis of gene expression data is the identification of groups of genes that manifest similar expression patterns. This translates to the algorithmic problem of clustering genes based on their expression patterns. Results: We presen

Tracing the Neural Carryover Effects of Interpersonal Anger on Resting-State fMRI in Men and Their Relation to Traumatic Stress Symptoms in a Subsample of Soldiers

Uncontrolled anger may lead to aggression and is common in various clinical conditions, including post traumatic stress disorder. Emotion regulation strategies may vary with some more adaptive and efficient than others in reducing angry feelings. However, such feelings tend to linger after anger provocation, extending the challenge of coping with anger beyond provocation. Task-independent resting-state (rs) fMRI may be a particularly useful paradigm to reveal neural processes of spontaneous recovery from a preceding negative emotional experience. We aimed to trace the carryover effects of anger on endogenous neural dynamics by applying a data-driven examination of changes in functional connectivity (FC) during rs-fMRI between before and after an interpersonal anger induction (N = 44 men). Anger was induced based on unfair monetary offers in a previously validated decision-making task. We calculated a common measure of global FC (gFC) which captures the level of FC between each region and all other regions in the brain, and examined which brain regions manifested changes in this measure following anger. We next examined the changes in all functional connections of each individuated brain region with all other brain regions to reveal which connections underlie the differences found in the gFC analysis of the previous step. We subsequently examined the relation of the identified neural modulations in the aftermath of anger with state- and trait- like measures associated with anger, including brain structure, and in a subsample of designated infantry soldiers (N = 21), with levels of traumatic stress symptoms (TSS) measured 1 year later following combat-training. The analysis pipeline revealed an increase in right amygdala gFC in the aftermath of anger and specifically with the right inferior frontal gyrus (IFG).We found that the increase in FC between the right amygdala and right IFG following anger was positively associated with smaller right IFG volume, higher trait-anger level and among soldiers with more TSS. Moreover, higher levels of right amygdala gFC at baseline predicted less reported anger during the subsequent anger provocation. The results suggest that increased amygdala-IFG connectivity following anger is associated with maladaptive recovery, and relates to long-term development of stress symptomatology in a subsample of soldiers